Complete mitochondrial genome of Morphostenophanes yunnanus (Zhou, 2020) (Insecta: Coleoptera: Tenebrionidae) and phylogenetic analysis

Abstract Morphostenophanes yunnanus (Zhou, 2020) is widely distributed from central to eastern Yunnan with distinct geographical variations in morphology. Beetles were collected in Manwan Town, and a mitochondrial genome sequence (GenBank accession number MZ298928) of this species was sequenced using the MGI-SEQ 2000 platform, assembled using NOVOPlasty v4.3.1, and characterized. The mitogenome was a circular DNA molecule of 15,690 bp with 64.710% AT content, which comprised 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region. The protein-coding genes showed the typical ATN (Met) and TTG (Met) start codons, except nad1 and cox1 (TTG as start codon), and were terminated by typical TAN stop codons. The maximum-likelihood polygenetic tree was generated using protein sequences of thirteen protein-coding regions of seventeen mitogenomes with mtREV + G + F + I with 1000 replicates under the Bayesian information criterion using MEGA 11, which showed that M. yunnanus was the most closely related to M. sinicus. This study provides essential genetic and molecular data for phylogenetic analyses of the genus Morphostenophanes.

Morphostenophanes yunnanus (Zhou, 2020) is widely distributed from central to eastern Yunnan with distinct geographical variations in morphology (Zhou 2020). It has gray-black shagreened mouth parts and dark-brown tarsi (Zhou 2020). M. yunnanus can be divided into seven geographical populations (Zhou 2020). The complete mitochondrial genome of M. yunnanus of the Manwan Town population was sequenced and characterized.
Adult M. yunnanus specimens were collected from Shuibatou Village (100.32963 N, 24.68407 E), Manwan Town, Yun County, Lincang City, Yunnan Province, China, on 24 September 2020, and deposited in the Insect Collection of Institute of Plant Protection, Guangxi Academy of Agricultural Sciences (http://www.gxaas.net/s.php/zwbhyjs/, Xuyuan Gao, gxy@gxaas.net) under the voucher number GIPP-20200924-001. Genomic DNA was isolated and subjected to paired-end sequencing (2 Â 150 bp) of 400 bp inserts using the MGI-SEQ 2000 platform. We produced $11.537 Gb of raw data, of which 11.301 Gb (97.76%) were clean high-quality data using SOAPnuke version 2.1.0 (Chen et al. 2018) with default parameters. The genome was assembled de novo using NOVOPlasty v4.3.1 (Dierckxsens et al. 2017) with default parameters and the mitogenome of Morphostenophanes sinicus (MW853764.1) (Bai, Gao, et al. 2021) as a seed sequence (no extending the seed directly). Three pairs of primers were designed corresponding to the cox1, nad5, and control regions (CR, also an AT-rich region) to verify the accuracy of the genome assembly. PCR products were sequenced from other M. yunnanus specimens using Sanger sequencing using an ABI 3730 automatic sequencer and assembled manually (Supplementary 1, 2), and only a few bases were found to be different (Supplementary 1).
For phylogenetic analyses, M. yunnanus mitochondrial PCGs and 16 other insect species were used to construct a maximum-likelihood (ML) phylogenetic tree. Protein sequences on thirteen protein-coding regions of their mitogenomes were aligned using MEGA 11 software (Tamura et al. 2021) with MUSCLE program (Edgar 2004) with default specifies. The ML Model with the lowest Bayesian Information Criterion (BIC) scores is considered to be the best. According to BIC (¼93859.54), General Reversible Mitochondrial (mtREV24) þ a discrete Gamma distribution (G) þ amino acid frequencies (F) þ evolutionarily invariable (I) with 1000 replicates was chosen to construct a phylogenetic tree using MEGA 11 (Tamura et al. 2021) (Figure 1). Sixteen mitogenomes of species belonging to the family Tenebrionidae and the mitogenome of Lepisma saccharina (Bai et al. 2020) were selected as the outgroup. The resulting polygenetic tree revealed that M. yunnanus was the most closely related to M. sinicus. Overall, our study provides insights into the mitogenome of M. yunnanus and essential genetic and molecular data for phylogenetic analyses of the species belonging to the genus Morphostenophanes.
Phylogenetic relationships were inferred using maximum likelihood and mtREV þ G (parameter ¼ 0.4735) þ F þ I (17.50% sites). The highest log-likelihood of the tree is À48283.84. The percentage of trees is shown below the branches in red. The complete mitogenome of M. yunnanus (GenBank accession number MZ298928) determined in this study is indicated in red.

Author contributions
Yu Bai analyzed the data, uploaded the analysis data, involved in certain tools for analysis, drafted the paper, and approved the final draft. Yang Kang performed the experiments. Lin Ye collected and analyzed data. Xuyuan Gao identified insects, contributed reagents/materials, involved in the conception and design of the work, performed the experiments, prepared figure, and approved and published the final draft. All authors agree to be accountable for all aspects of the work.

Disclosure statement
No potential conflict of interest was reported by the author(s).

Data availability statement
The genome sequence data that support the findings of this study are openly available in GenBank of NCBI at https://www.ncbi.nlm.nih.gov under the accession no. MZ298928. The associated BioProject, SRA, and Bio-Sample numbers are PRJNA729642, SRR14523735, and SAMN19159615, respectively.